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We present an integral formulation of the vacuum energy of Schrodinger operators on 
finite metric graphs. Local vertex matching conditions on the graph are classified accord- 
ing to the general scheme of Kostrykin and Schrader. While the vacuum energy of the 
graph can contain finite ambiguities the Casimir force on a bond with compactly sup- 
ported potential is well defined. The vacuum energy is determined from the zeta function 
of the graph Schrodinger operator which is derived from an appropriate secular equation 
via the argument principle. A quantum graph has an associated probabilistic classical 
dynamics which is generically both ergodic and mixing. The results therefore present 
an analytic formulation of the vacuum energy of this quasi-one-dimensional quantum 
system which is classically chaotic. 
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1. Introduction 

Quantum graphs have recently become popular as simple models of complex quan- 
tum systems as results derived on graphs have a good track record of providing 
insight in other areas. Some examples of fields where quantum graphs have been 
usefully employed include Anderson localization, photonic crystals, microelectron- 
ics, the theory of waveguides, quantum chaos, nanotechnology, and superconduc- 
tivity; see Refs. [1] and [2] for reviews. Here we consider a quantum graph model of 
vacuum energy. The vacuum energy of a star graph was first described by Fulling, 
Kaplan and Wilson where they note that for Neummann like matching conditions 
the Casimir effect is repulsive!^ Subsequently the vacuum energy of the Laplacian 
on a general metric graph was formulated as a periodic orbit su mil and in terms of 
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the vertex matching conditions by the authorsP^ One interesting feature of quan- 
tum graphs is that they possess a probabilistic classical dynamics that is chaotic. 
Such a mathematical model of vacuum energy for a spectrum of a chaotic system 
follows in the spirit of the Riemannium modei^ where the height up the critical line 
of the zeros of the Riemann zeta function provides the spectrum. Both the height 
of the zeros of the Riemann zeta function on the critical line and the spectrum of a 
generic quantum graph (with time reversal symmetry broken by a magnetic field) 
are distributed like the eigenvalues of a large random Hermitian matrix from the 
Gaussian Unitary Ensemble @ a typical feature of the spectrum of a system that is 
classically chaotic P 

In order to regularize the vacuum energy we employ the spectral zeta function, 

oo 
3=0 

where 7 e I is some spectral parameter and the sum runs over the spectrum 
of the Schrodinger operator < Eq ^ Ex ^ . . . ; the condition that the spec- 
trum be positive can be relaxed by introducing some additional technicalities.^ The 
zeta function is constructed from a secular equation whose positive roots kj satisfy 
Ej = fcjEI From the secular equation the zeta function is derived via the argument 
principle, as in Ref. [10] where the authors apply the zeta function to derive the 
spectral determinant of the Schrodinger operator on a graph with general vertex 
matching conditions. 

The vacuum energy is defined as 

£ c = lim^C(e- 1/2,0), (2) 

which is formally \ T^jLo t 1 1S an arbitrary parameter with the units of mass, 

which is included as £( s i7) wm generically have a pole at s = —1/2. Taking the 
limit e — > 0, 

S c = 1fPC(-1/2,0) + J f i+ln/1 2 ) Res C(-l/2,0) + 0(e) , (3) 



2 ,v ' ' ' 2 

where FP and Res denote the finite part and residue respectively. From the zeta 
function we obtain the following general result for the Casimir force on a bond j3 of 
the graph. 

Theorem 1.1. For the Schrodinger operator on a graph with local vertex matching 
conditions defined by a pair of matrices A and B, with AB^ = BA^ and rank(A, B) = 
2B, and potential functions Vb(xb) G C* 00 for b = 1, . . . , B the Casimir force on a 
bond j3 where suppl/g C (0,Lp) is given by 



f) 1 r°° f) 



a When the spectra are unfolded so they have mean spacing one. 
b The degeneracy of the roots and energy eigenvalues also coincide. 
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where F(it) = det (A + MM (it)) and J 7 ^ Dir is the Casimir force on the interval 
[0,Lp] with Dirichlet boundary conditions. 

M(it) is a 2B x 2B matrix denned in terms of the solution to a boundary value 
problem on the intervals [0, Lb]; see Section [3l The secular equation whose solutions 
k are square roots of the eigenvalues of the Schrodinger operator reads F(k) = 0. In 
the following we assume F(0) ^ 0, a condition that can be relaxed!^ The Casimir 
force on an interval with Dirichlet boundary conditions is defined in Eq. (|37|) . 

The article is organized as follows. Section [5] introduces the Schodinger opera- 
tor on a graph. Section [3] describes the secul ar e quation we use for the quantum 
graph. In Section 0] we recall our recent result for the zeta function of the graph 
Schrodinger operator. In Section[5]we evaluate the vacuum energy of the graph and 
prove Theorem ll.il 

2. Graph model 

A graph G is a collection of vertices v — 1, . . . , V and bonds b — 1, . . . , B, see Fig. 
[TJ Each bond connects a pair of vertices b = (v,w) and we use o(b) = v to denote 
the origin vertex of b and t(b) — w for the terminal vertex, b = (w, v) will denote 
the reversed bond with the origin and terminus exchanged. We consider undirected 
graphs with B undirected bonds, b and b will denote to the same physical bond 
with the label b referring to the bond with o(b) < t(b). m v denotes the degree or 
valency of v, the number of bonds meeting at v. To determine the valency b and b 
are a single bond. 




Fig. 1. A finite graph with 6 vertices and 7 bonds. 

In a metric graph each bond b corresponds to an interval [0, LJ. For the coor- 
dinate Xb € [0, Lb], Xb = at o(b) and Xb = Lb at t(b), Lb is the length of bond b. 
We also use the coordinate xj- = Lb — Xb which measures the distance to Xb from 
the terminal vertex of b. A function ip on G consists of a collection of functions 
{ipb(xb)}b=i,...,B on the set intervals. The choice of coordinate on a physical bond 
enforces the relation ipb(xb) = ^kbi x b) — ~ x b)- The Hilbert space of G is then 

B 

«=©i 2 ([0,id) • (5) 
f)=i 
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The Schrodinger equation on the interval [0,L&] reads, 

2 



h A b ipb{xb) + Vb(xb)ipb(xb) 

dxh 



k 2 M*b) ■ (6) 
= — At as the direction is 



A b is a vector potential on b and for consistency At = 

reversed when changing coordinated Matching conditions at the set of vertices are 
specified by a pair of 2B x 2B matrices A and B according to 



Aip + Bip = 



(7) 



where 



?P= (MO),---,iPb(0),MLi),--.,iPb(Lb)) , (8) 
ip= (0^(0),..., D Bl p B (0), Dj^x {L 1 ),...,D^I> B {L B )f . (9) 

Here, Db := g^- — iAb is the covariant derivative, so xp is the vector of inward 
pointing covariant derivatives at the ends of the intervals. The Schrod inge r operator 
is self-adjoint if and only if (A, B) has maximal rank and AB^ = BA^PJJ Such pairs 
of matrices classify the self-adjoint realizations of the Schrodinger operator. (See 
Ref. [12] for an alternative unique classification scheme for general vertex matching 
conditions.) 

For example, if we specify that the wavefunction is continuous at a vertex, 5-type 
matching conditions, a choice of matrices A and B that encodes for this is 



/-A 
-1 1 
-1 1 








and 



H = 



\ - 1/ 



Z 11 





Vo o 



1 x \ 


... o 



(10) 



0/ 



where A fixes the strength of the coupling at the vertex. If A = the incoming 
derivatives at the vertex sum to zero, Neumann like matching conditions. 

We further restrict the possible matching conditions to local matching conditions 
where the matrices A and B only relate values of functions and their derivatives on 
the intervals where they meet at a vertex and where the matrices are independent 
of the metric structure of the graph, namely the bond lengths and the spectral 
parameter k. Local matching conditions will be required subsequently to use the 
argument principle. 



3. Secular equation 

We first consider the simplest graph, an interval [0, L] with Dirichlet boundary 
conditions, i/j(0) — ip(L) = 0. Let f(x; —k 2 ) be a solution of the Schrodinger Eq. © 
such that /(0; — k 2 ) = 1 and f(L; —k 2 ) = 0. A second linearly independent solution 



C A vector potential Ab(xi,) that depends on xj, can be made constant by a gauge transformation. 
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is denoted f(x; —k 2 ) where x — L — x and / satisfies /(0; — k 2 ) = 1 and f(L; —k 2 ) = 
0. An alternative set of solutions are Uk(x; —k 2 ) and u(x; —k 2 ) satisfying it(0; — k 2 ) = 
and u'(0; — k 2 ) = 1. The function u is then, 

— / — ; 2\ f( X l—k 2 ) 

u{x >- k) = - r(L;-k 2 ) - 

k > is a square root of the eigenvalues of the Schrodinger operator iff 

u(L;-k 2 ) = -l/f(L;-k 2 ) = 0, (12) 

which defines a secular equation for the Dirichlet problem on the interval. For a 
graph with Dirichlet boundary conditions at the end of each interval [0, L b ] the set 
of intervals are effectively detached and the spectrum is the union of the spectrum of 
each of the intervals. Hence a secular equation for the graph with Dirichlet boundary 
conditions is 

B 

l[(-l/fl(L b ;-k 2 )) =0. (13) 

6=1 

If we now consider the Schrodinger equation on a general graph with matching 
conditions defined by the matrices A and B the component of a wavefunction with 
energy k 2 on bond b can be written 

Mxb, ~k 2 ) = c b f b (x b ; -k 2 ) e iA ^ + ^ f b {x t , —k 2 ) e iA ^ , (14) 

where f b is a solution of the eigenvalue problem on the interval [0, L b ] as defined 
previously. The vectors of values of the wavefunction and its inward covariant deriva- 
tives are then 

ip = (c!,...,Cb,Ci, ...,c w ) t and t/> = M{-k 2 )ip , (15) 

where M(—k 2 ) is the 2B x 2B matrix 

M ab = 6 a , b fi(0; -k 2 ) - 6 a l fi(L b ; -k 2 )e iAbLb . (16) 

Substituting in the vertex matching conditions we find that k 2 is an eigenvalue of 
the Schrodinger operator iff k is a solution of the secular equation 

F(k) := det(A + BAf(-fc 2 )) =0 . (17) 



Notice that the function F(z), where z = k + it, has poles along the real axis 
at points where l/f^(L b ; — z 2 ) vanishes, i.e. at the eigenvalues of the graph with 
Dirichlet boundary conditions at the ends of the intervals. 



4. Zeta function 

We construct the zeta function of the graph from the secular equation using the 
argument principle!^ For an interval with Dirichlet boundary conditions the secular 
equation reads u(L;—k 2 ) = 0. The function u(L;—k 2 ) has no poles and we can 
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represent the zeta function with a contour integral where the contour C encloses the 
positive real axis, namely 

Cuir(s, 7 ) = / (z 2 + l)- s ^-log (u(L- -z 2 )) dz . (18) 

Transforming the contour C to the imaginary axis z — it, Figure [2£b) 3 we obtain 

f°° l 

C Dir ( S , 7 ) = — / (i 2 -7)- s ^log( W (L;t 2 ))di, (19) 

a representation which is valid in the strip 1/2 < Jis < 1; see Ref . fTUl for details. For 
example, for a free particle where V(x) = 0, u(x;—k 2 ) = sin(kx)/k so u(L;t 2 ) — 
sinh(iL) ft. 




b) 



o zero of F(z) 
x pole of F( z) 



o zero of F(z) 
x pole of F( z) 



Fig. 2. The contours used to evaluate the graph zeta function, (a) before, and (b) after, the 
contour deformation. The two branch cuts are represented by hashed lines. 



For a graph with Dirichlet boundary conditions at the ends of all the intervals the 
zeta function is the sum of the Dirichlet zeta functions on the individual bonds. A 
graph with general local vertex matching conditions has secular equation F(z) = 
defined in Eq. (fT7|) . F has zeros on the positive real axis at square roots of eigenvalues 
of the operator and poles at eigenvalues of the graph with Dirichlet conditions. We 
assume we are in the generic case where these are distinct. The contour C is then 
chosen so that it encloses the zeros of F and excludes the poles on the positive real 
axis, see Fig. |2Ja). Transforming the contour to the imaginary axis, Fig. [U(b), we 
see that the zeta function has two components, 

B 

C(s,7)=Clm(s,7)+^CDir(s,7) ■ ( 2 0) 
b=l 

CDi r is the zeta function of the graph with Dirichlet boundary conditions at the ends 
of the intervals which comes from the poles of F and Cim is the contribution of the 
integral along the imaginary axis, 

Cim(s,7) = /V-7)- s -|logF(i*)di , (21) 

which converges in the strip 1/2 < 5Rs < 1. 
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In order to evaluate the vacuum energy we need to extend the range of validity to 
include s — —1/2. The restriction to 5fts > 1/2 comes from the asymptotic behavior 
of F as t — !• oo. The asymptotics of ]' h can be analyzed using the WKB methodp^l 
and we find 

oo 

ti(0;t 2 ) ~ -i + J>,,-(0)t-'' , (22) 

J'=l 

and f{,(Li,;t 2 ) vanishes exponentially. The functions Sb,j(xb) are determined by the 
recurrence relation 

Sb,j+i(xb) = - [ s' b Ax b ) +y^Sb,i(x b )sb,j-i(xb) ) , (23) 

where s b -i(x b ) = — 1, S6,o(^f>) = and Sb,i(x 6 ) = - V 6 (x b )/2. Using 
write the expansion of F(it) up to exponentially suppressed terms as 



F(it) ~ dot (A + M(-tl + D(t)) , (24) 

where 

oo 

D(t) = J2 * _i dia § {«i,i(0), • • • , WO), s Tj (0), . . . , «5j(0)} • (25) 
i=i 

This expansion of F(it) fits with the free particle case investigated previously in Refs. 
■5: and |6] where, in the absence of a potential, D(t) — and F(it) ~ det(A — iB). In 
general the expansion takes the form 



Let us denote by c/v the first nonzero coefficient in the expansion, so N = if 
Co = detB =/= 0, and by c^+j the second non-zero coefficient. Hence, 

logF(ii) ~ l 0g{cN t 2B - N ) + 2E±L + (t-^) . (27) 

i-)-oo C/\rt J 

Adding and subtracting these first two terms in the expansion of _F we obtain a 
representation of £im(s,7) convergent in the strip — (J + l)/2 < 5Rs < 1, namely 



sin7rs / >0 °^2 \-a d 
Cim s,7) = / (t 2 -7) t: 



log (^(ii)^- 2B / Cjv ) - £5tg 

V / CAT t J 



(It 



(2B~N)smir S _ s Jc N+J T(s + J/2) _ (2s+J) 

2™ 7 2 CiV r(s)r(j/2 + i) 7 ■ 1 j 

We must also represent the Dirichlet zeta function of each interval in a strip 
containing s = —1/2. The large t asymptotics of logit;,(i&; t 2 ) are given bjff^ 

logu 6 (i 6 ;t 2 ) - tL 6 - log(2t) + dfe t _1 + 0(t -2 ) , (29) 

t— >oo 
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where d b = | J^ b V(x b ) dx b . Subtracting and adding the first three terms we obtain 
the representation 



(t 2 - "f)- s ^ [log {u b {L b - t 2 )) - tL b + log(2t) - d b t- 1 ] dt 



T( S - 1/2) x , 1 _ s 2r(s + l/2) _x s 

2v7rr(s) 2 v 7rr ( s ) 



valid in the strip — 1 < 5Rs < 1. 



(30) 



5. Casimir effect 

Following Eq. @ , to evaluate the vacuum energy we identify 

,.\ ■ -2B _ C N+.J 

CN t J 



1 [°° d 

FPCim -1/2,0 = — / t- 

7T J 1 dt 



1 



f- 



dt 



log (F(it)) 

a pole i 
ResCim(-l/2,0) = 



dt . 



dt 



In the generic case when J = 1 there is a pole of ((s, 0) at s = —1/2 and 

CjV+l 
27TCJV 

For the interval [0,£f,] with Dirichlet boundary conditions 

I [°° d 



(31) 
(32) 



FP Cdu (-1/2, 0) = ~ / ^ [log («&(£&; t 2 )) - t£ 6 + log(2i) - d b r 1 ] dt 

-\J ' *^ [log (u b (L b ;t 2 ))] dt , (33) 

and Res£j3 ir (— 1/2, 0) = d b /n. Combining these results and using Eqs. ^ and (|2U)) 
we have obtained a representation of the vacuum energy of the graph Schrodinger 
operator. 

If we assume that the potential on a bond is compactly supported and consider 
the Casimir force on that bond we can obtain a more explicit description of the 
Casimir effect on the graph. For a bond j3 let 

supply C (0,Lp) , (34) 

then the Casimir force on the bond j3 is finite and can be evaluated without am- 
biguity as the asymptotic behavior of F(it) is independent of Lp. We define the 
Casimir force on bond j3 for the compactly supported potential as F@ — —-t^£ c , 
then 



7T/3 _ tt/3 



1 

27 



dL, 



\ogF(it) dt 



Dir ^ S ^ ne f°rce on the interval [0, Lp] with Dirichlet boundary conditions, 



i 

'2tt 



dL 



(35) 



(36) 
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Hence we obtain Theorem ll.il 

For example, with no potential on the graph 



r/3 1 r 9 ,_Y sinh 0M) 



^v=-^x ^(-^j*— (37) 



a result that can also be obtained directly from the spectrum Ej = Vj-j of 



_>2 

J 3 ~ V L 

the Laplace operator on an interval with Dirichlet boundary conditions, where 

c Dir M) = (£r 2s CR(2, S ). 

For a free particle we can also write f b explicitly, 

2 sinfc(i 6 - z 6 ) 2 -tcosht(L 6 - x b ) 

f b (x b ;-k ) = — — , hence f b (x b :t ) = — — . (38) 

sin kL b sinn i.Lj, 

The matrix M then has the form 

M(t a )=tf-™Sf CSC ^ T ) , (39) 
v ; V cschiL -cothiL / v y 

where cothiL = diag{coth£Li, . . . , coth t Lb} and cschiL is defined similarly. If 
we choose a particular graph like the star graph, Fig. |3l with Dirichlet boundary 
conditions at the vertices of degree one and <5-type coupling at the central vertex, 
such matching conditions can be encoded in the 2B x 2B matrices, 

A s J \0 B 5 

with As and Ms defined as in Eq. (|10p. Substituting in (|17l) and evaluating the 
determinant we find 

B 

F(it) = -A-^coth£L 6 , (41) 

b=l 



(40) 



and the Casimir force on a bond f3 is given by, 



i />oo a 




Fig. 3. j4 star graph with 5 bonds. 

If the condition that the potential Vp is compactly supported is relaxed, finite 
renormalization ambiguities in the Casimir force remain. These are the standard 
ambiguities as they o ccu r in self-interacting (A$ 4 )-theories and they can be dealt 
with in a routine wayP^ 
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